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, Abstract 

<^ We determine the topological susceptibility in the SU(3) pure gauge theory. We perform 

a series of high-statistics lattice studies and take the combined continuum and infinite 
volume limit. We find Xtoprf = 0.0524(7) (6), which translates into xllp = 193(l)(8)MeV 
with the second error exclusively due to the intrinsic scale ambiguity. 
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^ '. 1 Introduction 
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In QCD with N c colors and Nf light dynamical quarks there are two notions of the topological 
susceptibility, defined as the second moment of the global topological charge distribution. 
On the one hand, the actual topological susceptibility Xtop° shows a clear sensitivity on the 
Ph! dynamical (sea) quark masses [Tj. This property renders it an ideal vacuum diagnostics tool, 
as emphasized in [2] and exploited in a number of recent studies [31 HJ G2, El E El 13 [10] • 

On the other hand, the quenched topological susceptibility x™ is t ne quantity which is 
linked via the famous Witten-Veneziano formula (TTj [T2] 



X 



X™ = W (M 2 V ,+M 2 V -2M 2 K ) 



to the excess of the 7/ mass over the pseudoscalar octet masses, with a proportionality factor 
which contains the pseudoscalar decay constant^ in the chiral limit. This relation is supposed to 
hold at the leading order in an expansion in 1 / N c (for a different viewpoint see [131 El 03] ) • To 
be precise, this quenched susceptibility is the susceptibility of the underlying S*£7 (N c ) Yang-Mills 
(YM) theory, and this is what makes ([1]) appealing from a theorist's viewpoint - it relates two 
different theories. Specifically, on the lattice one can measure the l.h.s. of ([1]) for several N c and 
evaluate the r.h.s. for various (N c , Nf) combinations, and finally check whether the agreement 
is parametrically controlled by 1/N C or Nf/N c . 

In this paper we elaborate on the first step in this program - we determine with unprece- 
dented precision the quenched topological susceptibility for the case A^ = 3. We begin with an 



use the Bern normalization where FP hys = 92.4(3) MeV and F = 86.2(5) MeV in the chiral limit. 
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exposition of the main avenues towards defining a topological charge on the lattice and how 
one extracts the topological susceptibility. The next two sections contain details of our lattice 
simulations and of our combined continuum and infinite volume extrapolation. Having a result 
in terms of the Sommer radius r [15], the latter needs to be identified with a length-scale 
in fm, and we discuss both the input that goes into such an identification and the remaining 
ambiguity. In the concluding section we compare our result to other recent determinations of 
Xtop = Xtop 1 - Details of a new parameterization of r (/3) have been arranged in an appendix. 



2 Topological charge definition 

In the continuum the topological charge of a given gauge background is defined as 
q = J dx tr(F fllJ (x)F fllJ (x)) = J dx tr(e IJ , l/ap F fllJ (x)F ap (x)) 



where F fJiU = F^X a /2 is the field strength tensor. For toroidal space-time geometry (|2j) is integer 
and linked to the index of the Dirac operator D via the Atiyah-Singer theorem [17] 

q = n + — n_ (3) 

where n± denotes the number of zero modes of D with positive or negative chirality. 

On the lattice the definition of the topological charge is not unique. Aiming for the glu- 
onic side of ([3]), one may choose any discretization of F^ which has the correct perturbative 
continuum limit, and form the so-called "naive" (or unrenormalized field-theoretic) charge 

9nai[£/] = —^2 E tr(e^ CTp F^(x)F CTp (x)) = — £ tr(F 12 F 34 + F 13 F 42 + F 14 F 23 ) . (4) 



In general this definition does not lead to an integecl charge. 

On the other hand, the fermionic side of ([3]) always yields an integer answer. This is most 
straightforward with the massless overlap^ operator D = p[l + Z^ w ,-p(-Dw,-p^ ) w-p) _1 ^ 2 ] [22] 
or with any other Dirac operator D which satisfies the Ginsparg- Wilson relation [23]. In this 
case the index can be written in the closed form [211 [25] 



qiorP] = -^ tr (75^) • (5) 

With a non-chiral Dirac operator (e.g. -Dw) explicit mode counting prescriptions may be set 
up. Below, the only point which matters is that (7f er [£7] is necessarily an integer. 
The continuum topological susceptibility at zero virtuality is defined as 

Xtop = Jim %^ (6) 

and this shows that a finite volume is mandatory for the definition. Still, contributions to (q 2 ) 
which grow less than linearly in the 4D volume V create finite volume effects in Xtop- 



2 There are two known exceptions, the Liischer [18] and the Phillips-Stone [19] definitions of qhh[U}. 
3 It does not matter whether one uses the Wilson operator Z?w or another doubler-free kernel [2Ql [21] 
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On the lattice the details of the topological charge definition reflect themselves in the precise 
form of the latticized version of (El) [26]. In case one starts with the gluonic definition (jlj), the 
traditional approach has been to form the bare susceptibility (q^/V which is then subject to 
both additive and multiplicative renormalization [27] 128 ] 29|. On the other hand, starting from 
the fermionic definition (jSJ), one just forms (qf er )/V, since is already a renormalized charge. 
Likewise, if we first compute a renormalized (integer) field-theoretic charge, the susceptibility 
based on it will not require any further renormalizatiorfl 

In this paper we investigate the utility of such a renormalized gluonic charge definition for 
a precision measurement of the topological susceptibility in the pure SU (3) gauge theory. We 
start from the standard "clover- leaf" definition of F fiv (x) (it uses the average of the antiher- 
mitean part of the plaquette U^ u in x,x—p,,x—u,x—p, — u) based on HYP smeared [M] gauge 
links. Plugging this into (jlj) we have the bare charge qw^] which is a real number. One of our 
g nai distributions (/? = 6.0, 12 4 , 3 HYP steps) is shown in the left panel of Fig.[TJ Thanks to CP 
symmetry only a multiplicative renormalization applies, and we opt for a non-perturbatively 
defined Z-factor. Still, there are various possibilities, and we choose a strategy which makes 
use of the fact that on fine enough lattices the overall distribution of q na i tends to cluster near 
integer values (cf. Fig. [I]). We find Z as the solution^ of 

min(x 2 ) where \ 2 = Yl ( Z( lnai[ U ] - round (Zg nai [[/])) (7) 
z>i u 

and use it to define, via rounding to the nearest integer, the renormalized field-theoretic charge 

g ren [^] = round (Zg nai [[/]) (8) 



4 Note that the round-to-integer operation ((HJ) brings in a global element, and our q rcn is not given by the 
integral over a local charge density. Accordingly, there is no conflict with the result by Stamatescu and Seiler 
that in general q(x) mixes with the identity and the correlator J q(x)q(0) dx has a contact term 30 . Any overlap 
based charge definition sticks out in the sense that it is an integral of a local density and avoids the mixing with 
the identity [31j [32l [33] . The argument for the absence of additive mass renormalization is then standard - in 
the zero-charge sector the susceptibility vanishes exactly. Note that it is essential to define the sectors with the 
same operator that is used in the %top estimator - if one uses one overlap charge to define the sector and another 
overlap charge to measure (q 2 ), then even the overlap based susceptibility has an additive renormalization. 

5 The restriction Z > 1 is a technical aspect of the minimization procedure to avoid the global minimum 
\ 2 = at Z = [3S] . In perturbation theory one finds Z = 1 + const g 2 with const > ■ 
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which, by construction, is an integer. This charge definition has already been used in [35] . 

The alert reader might be surprised by our frequent use of expressions like "we choose" or 
"we opt for" in this passage. Indeed, there is a huge amount of freedom in how one attributes an 
integer charge to a lattice configuration. However, according to the standard scaling hypothesis 
by Symanzik this ambiguity reflects itself in different 0(a 2 ) cut-off effects of observables built 
from the topological charge. It has been checked that the "disagreement rate" between any 
pair of charge definitions quickly vanishes with (3 — > oo [37J and that the gluonic or fermionic 
charge definition yields the same continuum limit for the topological susceptibility [29JI3EI 39J. 
Clearly, this is not a mathematical proof, but it is worth emphasizing that this is exactly the 
behavior that one expects to see, if one is in the Symanzik scaling regime. 



3 Lattice simulations 

Our goal is to perform a series of simulations in a fixed physical volume that will allow us to 
determine the topological charge distribution in the continuum limit (in that volume). This 
will be complemented by a second series of simulations (at a fixed lattice spacing) to assess 
possible finite volume effects. 

We use the MILC code [ID] to produce the SU(3) gauge ensembles. We choose the Wilson 
gauge action and run parameters as detailed in Tab. [I] (for the scaling series) and Tab. [2] (for 
the volume extrapolation series). The scale is set via the Sommer radius r [IS], based on a 
new parameterization of data from [H] as specified in the appendix. The physical value is a 
separate topic that will be discussed along with the presentation of the final result. 

We start with a closer look at Tab. [I] Throughout, our runs are designed to yield ~ 10 5 
measurements. Anticipating that the autocorrelation time of the topological charge rapidly 
grows with /3, we have increased the number of updates, n sepa , between adjacent measurements. 
Still, it turns out that the measured integrated autocorrelation times are somewhat larger on 
the finer lattices. The renormalization factor Z for g na i with 3 HYP steps, as defined in ([7]), 
seems to go monotonically towards 1 with increasing /3, as expected. Based on it we determine 
<l = ( lren defined in (JS]). As a first check, we measure |(<?)|, and this moment is consistent with 
zero on all lattices. The observable of interest, (q 2 ), is measured with 1% statistical accuracy (or 
better), throughout. We have also checked that using a 2 HYP or 4 HYP charge would change 
the suscpetibility by an amount which is an order of magnitude smaller than the statistical 
error of the (q 2 ) given. Our data are precise enough to evaluate the fourth moment of the 
distribution, and the kurtosis (q A )/(q 2 ) 2 — 3 turns out to deviate from zero, for most lattices 
in the scaling series, by about 10a. It seems independent of the lattice spacing and it is thus 
natural to ask whether a non-vanishing (g 4 ) / (q 2 ) 2 — 3 is a finite volume effect. 

This brings us to discuss the second series of runs, as detailed in Tab.[2j Potential finite 
volume effects are commonly attributed to infrared physics; we thus restrict ourselves to a single 
coupling, (3 = 6.0. Again, we aim for O(10 5 ) measurements per run. For larger volumes the 
"crest and valley" structure in Fig.[T]becomes less pronounced^, and we pursue the analysis with 
the Z-factor determined on the (6.0, 12 4 ) lattices. Now one expects the autocorrelation times to 
be independent of the volume. This happens to be true, except in the L/a = 10 case, where we 
find severe finite volume effects in (q 2 ), too (see below). Again is basically consistent with 



6 Here it is understood that the total statistics is kept fixed, so that the number of configurations in a given 
sector drops. Via appropriately increasing the statistics, Z can be determined in a volume as large as desired. 
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zero. The second moment (q 2 ) is measured, as before, with 1% statistical accuracy (or better) 
on all lattices. The main point of this series is clear evidence that the kurtosis (<? 4 )/ {(f) 2 — 3, 
which quantifies the deviation from a Gaussian distribution, tends to zero with V ^oo. 

This observation lets us plot the charge histograms for two lattices from the scaling series 
in Fig. [21 It turns out that the relative weight of the toplogically non-trivial sectors is almost 
consistent with a Gaussian form, while the q = sector shows a clear excess. It is thus tempting 
to define a new (g 2 ) g ^ as the width of this Gaussian, where the fit takes only the nontrivial 
sectors into account. In fact, due to CP symmetry the distribution is even, and we produce a 
histogram of \q\. The \q\ > 1 sectors are then fitted to a half-Gaussian. This procedure stabilizes 



L/a 



5.8980 
10 



6.0000 
12 



6.0938 
14 



6.1802 
16 



6.2602 
18 



6.3344 
20 



^sepa 
^-conf 

z 

nm(q) 

r int (sign(g)) 

<<zW> 2 -3 
(<?W>-3(g 2 > 
(<? 4 >-3(g 2 ) 2 



(A 



10 
100010 
1.2902 
0.93 
0.61 
0.89 
0.004(5) 
1.695(9) 
0.214(21) 
0.363(36) 
0.616(61) 
1.702(12) 



10 
101600 
1.2353 
2.18 
1.09 
2.05 
0.007(8) 
1.592(11) 
0.238(24) 
0.378(38) 
0.602(62) 
1.610(13) 



10 
103705 
1.1993 
4.83 
2.44 
4.42 
0.007(11) 
1.490(15) 
0.244(27) 
0.363(40) 
0.541(60) 
1.526(22) 



50 
101710 
1.1742 
2.27 
1.19 
2.10 
0.007(8) 
1.465(10) 
0.227(21) 
0.333(32) 
0.487(46) 
1.486(13) 



50 
112222 
1.1570 
4.65 
2.45 
4.25 
0.021(10) 
1.427(14) 
0.241(30) 
0.344(43) 
0.490(63) 
1.461(21) 



100 
105314 
1.1430 

4.61 

2.52 

4.27 
0.010(10) 
1.381(14) 
0.204(32) 
0.282(44) 
0.389(61) 
1.400(19) 



Table 1: Summary of the scaling series of runs. Unless stated otherwise, q is the renormalized 
(i.e. integer) gluonic topological charge g ren with 3 HYP steps. The subscript u q ^ 0" means 
that the distribution of \q\ with q^0 has been fitted to a half-Gaussian. Statistical errors have 
been estimated via a jackknife with the blocklength set to round (10ri nt (g)). 
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0.238(24) 
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5.202(51) 


8.242(54) 



Table 2: Summary of the volume extrapolation series of runs. Throughout, the charge renor- 
malization factors of the (6.0, 12 4 ) lattices have been used. For details see caption of Tab. 1. 
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Q 
P 

L/a 


o.UUUU o.Uyoo o.ioUz o.zoUz o.oo44 
12 14 16 18 20 


\Q )\ r o/L) 

(?Vo(r /L) 4 


oooU^oUJ oi)AO[il) ooZo^o I ) uvtoyvt) 04yU^ooj 
6401(69) 6066(95) 5907(66) 5808(92) 5568(84) 



Table 3: Summary of the data from the scaling series of runs, in units of 10 5 . Here, an extra 
0.68% error has been added w.r.t. what is reported in Tab. 1. 



L/a (13 = 6.0) 


10 12 14 16 18 


(<f)(r /L) 4 
<<?Vo(V£) 4 


4834(42) 6330(42) 6464(43) 6540(43) 6468(42) 
6428(79) 6401(53) 6458(56) 6545(64) 6474(42) 



Table 4: Summary of the data from the infinite volume extrapolation series, in units of 10 5 . 
Here, no extra error has been added w.r.t. what is reported in Tab. 2. 



the contribution from the tails (for this reason a similar "improved estimator" was used in |42j). 
but it also introduces a model dependence. 



4 Continuum and infinite volume extrapolations 

For the continuum and infinite volume extrapolation of the quenched topological susceptibility 
we convert our data to physical units by multiplying them with the appropriate power of r /a. 
In the scaling series the uncertainty of L/r reflects itself independently for each datapoint. 
Assuming that L/r is known (at each (3) with 0.17% precision (see Appendix) an extra 0.68% 
error needs to be included in Tab. [31 By contrast, in the volume extrapolation series the 
uncertainty of L/r is correlated, and no extra error has been included in Tab. HI 

For the continuum extrapolation in a fixed physical volume we plot our data against (a/r ) 2 . 
This power of the lattice spacing is motivated by a Symanzik analysis (both the Wilson action 
and our topological charge operator contain dimension 6, but no dimension 5 operators). In 
principle, the uncertainty of r /a also leads to a horizontal error bar, but its effect is negligible, 




Figure 2: Distribution of q rcn for the (6.0, 12 4 ) [left] and (6.3344, 20 4 ) lattices [right]. The fit 
to a Gaussian form excludes the q = sector; the excess gives rise to the kurtosis (<? 4 )/(<? 2 ) 2 — 3 
reported in Tab. 1, an effect that goes away in the infinite volume limit, as evident from Tab. 2. 
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Figure 3: Continuum extrapolation of Xtop^o m a fixed physical volume, V = (2.2394r ) 4 . 
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Figure 4: Infinite volume extrapolation of Xto P r o> versus r o/L [cf. and (r /L) 4 [cf. (1TU1) ] 
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Figure 5: Continuum extrapolation of (<? 4 )/ (q 2 ) 2 — 3 in a fixed physical volume, V = (2.2394ro) 4 . 
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Figure 6: Left: Infinite volume extrapolation of the kurtosis (q 4 ) / (q 2 ) 2 — 3 in log-log form; this 
quantity vanishes with V — > oo. Right: The reduced moment ratio (g 4 ) / (q 2 ) — S(q 2 ) versus 
r A /V. This plot includes constant fits to the last 2, 3 or 4 datapoints and a linear fit to all 5 
points. We cannot discriminate a finite from a vanishing limit of this quantity with V — >oc. 
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and we shall omit it. The resulting fits are shown in Fig.0 The observables (q 2 )(ro/ L) 4 and 
(<? 2 ) <?^o( r o / L) 4 extrapolate in a parallel manner, giving the continuum values 0.05092(71) and 
0.05205(71), respectively. The difference, if evaluated inside a jackknife, is significant. 

This brings us to the discussion of finite volume effects. Standard reasoning suggests that 
the dominant finite volume corrections in the YM theory come from glueball states travelling 
around the box. Since euclidean Green's functions fall off exponentially, this would imply 

Z V (L) = Z u (oo) (l + const e~ MaL + ...) (9) 



in close analogy to relative finite volume effects in QCD being proportional to exp(— M n L). 
Simple models for the YM vacuum (random population of unit cells by instantons or anti- 
instantons) suggest no finite volume corrections for the susceptibility and a 1/V type kurtosis^. 
As a result, one may consider the alternative form of the partition function 

Z V (L) = Z v {po) (l + const/L 4 + ...) . (10) 



The data of Tab.H] are shown in Fig.U If we assume exponential finite-volume effects, the 
data for the first estimator may be fitted all the way out to L/a = 10. The difference of the 
fit curve at L/a = 12 and at L/a = oo suggests that the continuum result for (q 2 )(r / L) 4 is 
suppressed by the amount 0.00175, and this gives a corrected {q 2 )(ro/L) 4 = 0. 05267(71). With 
the other estimator the fitted value for "const" in (jUJ) is consistent with zero, and we stay with 
(q 2 ) q ^o( r o/ L) 4 = 0.05205(71), as quoted above. If we assume 1/V- type effects, finite-volume 
corrections are hard to quantify (at least the volume-scaling regime with the first estimator is 
smaller). Therefore, we decide to stay with the extrapolation via (jUJ). Taking the average of 
both estimators but keeping the full difference as a theoretical uncertainty, we consider 

Xto P r 4 = 0.05236(71)(62) (11) 

our final result, where the first error is statistical and the second systematic. 

Our data for the kurtosis (q A )/(q 2 ) 2 — 3 scale perfectly (Fig.EJ), and the infinite volume 
extrapolation is consistent with both an exponential and a power law behavior. We find that in 
the YM theory a non-zero kurtosis is a pure finite volume effect and has nothing to do with the 
discretization (Fig.EJ). Since our "independent cell" model of footnote 7 and large N c arguments 
suggest that (g 4 ) / '(q 2 ) — 3(q 2 ) has a finite limit with V— >oo, we plot this quantity versus 1/V 
in Fig. El too. It turns out that our data are consistent with the value 0.276(84) given in |43j, 
but we cannot rule out the possibility that the infinite volume limit might actually vanish. 



5 Conversion to physical units 

Having an unambiguous result for the topological susceptibility in the combined continuum 
and infinite volume limit in units of the Sommer scale r^ 1 [16J, we are left with the question 
which value [in MeV] the latter should be identified with. 

In QCD r is a well-defined quantity. In other words, one could measure aM p , aM n and 
aMx in a 2 + 1 flavor simulations in which their ratios are always adjusted to the respective 

7 With independent fluctuations in unit cells V\, Vn this is generic: (v 2 ) = ((v\ + ... + vn) 2 ) = N (v 2 ) and 
(^ 4 ) = Hitf) + 6 E^-( I/ f)( J/ |) = N(vf) + 3N{N - l){v 2 )(P2), where we use that all odd moments vanish. 
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Xt£ [MeV] xt'ot [MeV] 



Ref. [49J (2001) 
Ref. [43] (2002) 
Ref. [39] (2003) 
Ref. [50] (2003) 
Ref. [42] (2004) 



0.072(7) 

0.057(3) 

0.055(10) 

0.059(5) 

0.059(3) 

0.0524(9) 



209(5)(8) 



198(3)(8) 191(5)(8) 
193(1)(8) 



197(3)(8) 180(2)(8) 
195(9)(8) 188(12)(8) 
198(4)(8) 



this work 



Table 5: Summary of the quenched topological susceptibility (in the continuum) as determined 
in some recent studies. The results Xto P /o" 2 = 0.0355(33) [49] and Xto P /o" 2 = 0.0282(12) [43] have 
been converted to r^ 4 units by means of cr 1 / 2 ^ = 1.193(10) [5T] . In the second column the 
conversion to MeV has been done with our choice ro = 0.49fm, while the third one contains the 
values given in the respective papers, but with our standard scale setting error throughout. 



experimental ratios. By considering the continuum limit of r M p , one would have the correct 
physical value of r , and there is nothing wrong with using this value also in pure YM theory. 
The original estimate r = 0.5 fm [16J for the outcome of this procedure has been superseded, 
more recently, by values like r = 0.467(6) fm [44J|45]. Alternatively, one could set the scale in 
the quenched theory via r fx = 0.4146(94) [46], using fx = 160(2) MeV from experiment. This 
has the same effect as attributing r = 0.512(12) fm. To encompass this spread we set the scale 
via r = 0.49fm and add a 4% error to reflect the intrinsic scale ambiguity in pure YM theory. 
With this choice for r our value Xtop^o, as quoted in the previous section, amounts to 



where the first error bar contains all statistical and systematic uncertainties of our calculation, 
and the second one reflects the scale setting ambiguity in a theory which is not full QCD. 

6 Summary 

In this paper we have performed a precision study of the topological susceptibility in pure SU (3) 
Yang Mills theory in the combined continuum and infinite volume limit using a field-theoretic 
definition of the topological charge. 

Comparing our final result to other recentjfl determinations [49], E3] [39], [50l |42] (see Tab. [5] 
for details) we see a downward trend over time. Our result in r$ 4 units is substantially more 
precise, since we work much closer to the continuum (our /3 = 6.0 lattices are the coarsest ones 
used in the continuum extrapolation) and due to the large statistical The conversion to MeV 
is, of course, limited by the basic scale setting ambiguity in the quenched theory. In any case 
our final result ({12]) supports the Witten Veneziano scenario for the origin of the rj mass. 

Acknowledgments: We thank Ferenc Niedermayer for useful discussions. Computations 
were carried out on ALiCEnext and a Pentium4 cluster at the University of Wuppertal. This 
work was supported by the Swiss NSF and under the Hungarian grant OTKA-AT049652. 

8 We like to draw the readers attention to two studies [47] 08] with a fermionic charge where no continuum 
limit has been taken. Note that the charges Qi in [49] and Q g in Fig. 5 of [39] are similar in spirit to our q len . 

9 There is a similar high-statistics study, based on an overlap charge [52] . They have no continuum limit yet 
for Xtop, but they see a deviation from zero in (g 4 ) — 3(q 2 ) 2 or the kurtosis in a fixed volume. 




193(1)(8) MeV 



(12) 
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Figure 7: Old polynomial x 2 /6 = 1.39 and new rational x 2 /6 = 0.91 fit to the data of Ref. [UJ. 
Asymptotic freedom demands that the curve eventually becomes linear with slope 47r 2 /33, and 
this constraint is built into the new rational ansatz (fl4|) . 



Appendix 

In Ref. jH] Necco and Sommer use a polynomial approximation for log(ro/a). Specifically 

log(r /a)| Ns = c + c x {(5 - 6) + c 2 (/3 - 6) 2 + c 3 (/3 - 6) 3 (13) 

with Co = 1.6804, c\ = 1.7331, c 2 = —0.7849, C3 = 0.4428 is suggested in order to generate matched 
lattices for 5.7 < (3 < 6.92, i.e. combinations of ((3,L/a) with fixed L/r . We have repeated 
the analysis with their data [UJ and find essentially the same Co.. .3, with % 2 /d.o.f. = 1.39. 

Here, we wish to explore a rational approximation of log(r /a) which has the merit of being 
consistent with perturbation theory. In the weak coupling regime a oc exp(—l/[2j3 g 2 }) with 
the universal coefficient Po = l/ (4n) 2 ■ [llA^/3 — 2Nf/3]. In SU(3) gluodynamics one has thus 
acxexp(— 167r 2 /132 • (3) or log(r /a) = 4-7r 2 /33 •/? + ... and we employ the ansatz 

log r /a = — ■ 14 

new 33 dzl p + 04/ p l 

where this constraint is built in. Using again their data we find that the best fit is given through 
di = -8.2384,d 2 = 15.310, d 3 = -2.7395,d 4 = -H.526 with x 2 /d.o.f. = 0.91. 

It is straightforward to find the /^-values which match a simulation with (/?, L/a) = (6.0, 12). 
The results are given in Tab. El both for ([TBI) and for (fill) . It turns out that the matched (3- 
values are more or less the same with either formula and the pertinent L/tq (and hence our 
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L/a 


10 12 14 16 18 20 


a 
P 

L/r 


o.oyyo o.UUUU o.Uyzo o.l/yU o.zoUl o.oooz 
2.2357 2.2356 2.2357 2.2357 2.2356 2.2357 


a 
P 

L/r 


o.syou o.uuuu o.uyoo o.iouz o.zouz 0.0044 
2.2394 2.2394 2.2394 2.2395 2.2393 2.2393 



Table 6: Upper part: (/?, L/a) combinations which match (6.0, 12) [via achieving L/r = 2.2356 
as accurately as possible], based on the interpolation formula fll3p of [JT]. Lower part: Same 
[via achieving L/tq = 2.2394 as accurately as possible], based on the interpolation formula f|T4|) . 



estimates for the physical box-length) differ by 0.17% only. This should not come as a surprise, 
since the two curves in Fig. [7] are rather close. Still, since the new interpolation has a better 
theoretical foundation and the xVd.o.f. is lower, we choose the new parameterization (jl"4"l) . 
Nonetheless, it is clear that the original parametrization f|T3|) remains a legitimate and good 
one in the interval in which it is given. 
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